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CHAPTER I 


INTRODUCTION 

Nickel base single crystal superalloys, because of their superior 
high- temperature properties, have attracted considerable interest for 
use in gas turbine engines. The initial development of these alloys as 
turbine blade and vane materials began in the 1960’s at Pratt and 
Whitney Aircraft with the alloy PWA 1 409 which is the single crystal 
version of Mar-M200, [1]-[5]. Since that time a variety of similar 
alloys have been developed, both by Pratt and Whitney and by other 
aircraft engine companies. These include the Pratt and Whitney alloy 
PWA 1480, the AiResearch alloy Mar-M247 and the General Electric alloy 
Rene-N4, which is the alloy of interest in this research. 

During the manufacture of single crystal turbine blades the [001 ] 
crystal orientation is the natural growth direction. The [001] axis is 
parallel to the span of the blade which i3 also in the direction of the 
centrifugal loading. At this orientation considerable improvement in 
creep and fatigue resistance is obtained over conventionally cast 
turbine blades. There is some variation in the [001] orientation from 
one blade to another due to the manufacturing process. The transverse 
orientation of the crystal is generally not controlled, but it can be 
controlled by using a seed crystal. At the present time it is not 
possible to adequately evaluate the benefits of controlling the 
transverse orientation. 

One of the basic needs is the development of better mechanical 
analysis methods for use in the design of single crystal engine 
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components. Specifically, it is necessary to develop a model to predict 
the inelastic constitutive behavior of the material and to incorporate 
the model into a general purpose finite element code. The purpose of 
this research is to satisfy this need. Presented in this study is a 

constitutive model for Rene N4 at 760°C and the implementation of this 
model in a nonlinear three dimensional finite element code. 

1 .1 Alloy Chemistry and Structure 

Rene N4 was developed by General Electric specifically for use as 
a turbine blade and vane alloy. The nominal chemical compositions of 
Rene N4 and PWA 1480 are presented in Table I, along with the two 
specific compositions examined in this study. References [6] and [68]. 
The elements C,B, and Zr which are typically included in polycrystalline 
high temperature, nickel base superalloys for* the purpose of 
strengthening the grain boundaries are omitted. 

All of the modern single crystal materials are two phase alloys 
with a large volume fraction of Y* phase. The Y' precipitates have L1 2 

type crystal structure and are interspersed in a coherent face centered 
cubic Y solid solution, Figure 1. The strength of the alloy is a 
function of the Y’ size and the percentage of Y', [7], [15]. 
Experimentally it has been determined that the peak creep strength is 
achieved with a volume fraction of Y 1 of about 60J. In high volume 
fraction Y • alloys deformation must occur by shearing of the 
precipitates. Much of the behavior of the alloy can be explained on 
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this basis. As a result, considerable research has been devoted to 


crystals with LI 2 structure, [8] - [14]. 


TABLE 1 . 

Composition and Heat Treatment of Rene N4 and PWA 1480 



Element, WT. % 

A1 

Ti 

Ta 

Nb 

Cr 

Mo 

TI 

Co 

C 

B 

Zr 


^ENE N4 
NOMINAL 

o 

m 

m 

D 

9.25 

1.5 

a 

sa 


■ 

- 

C 

RENE N4 
COMPOSITION A 

3.77 

4.24 

3.96 

B 

9.26 

1.60 

5.88 

7.53 

Bf 

fl 

- 

d 

RENE N4 
COMPOSITION B 

3.6 

B 

4.1 

B 

8.7 

1.6 

B 

B 

■ 

■ 

- 

b PWA-1 480 

B 

1.5 

12 

B 

10 

fl 

g 

B 

- 

■ 

■ 


a Heat Treatment: 1 260°C/2h/ gas quench ♦ 1 080°C/4h/ air cool + 
900°C/1 6h/alr cool. 

b Heat Treatment: 1 288°C/4h/air cool ♦ 1079°C/4h ♦ 871°C/32h. 

Q 

Actual composition for specimens reported in References [6], [37]. 

d Actual composition for specimens reported in Reference [68] and 
designated In VF 317. 
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b) Gamma Prim Morphology at High Magnification (10,000X) 


FIGURE 1. MICROSTRUCTURE OF FULLY HEAT-TREATED RENE 
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1.2 Mechanical Response of the Alio' 


Elastic strains are associated with crystal lattice distortions 
and are fully recoverable whereas inelastic strains result from the 
movement of atoms within the lattice and are not recoverable. The 
effects of Inelastic strains on the elastic properties are assumed to be 
negligible. 

Rene N4, as well as other single crystal alloys, exhibits cubic 
symmetry in the elastic range. Reference [78]. There are three 
independent elastic constants and in the principal material directions 
the elastic constitutive equation is written as 



The components of the stress and elastic strain tensors are denoted by 
and e^. Since the elastic modulus, E, shear modulus, G, and 

Poisson’s ratio, v, are independent elastic constants the terms of the 
compliance matrix depend on orientation. Rotation of the elastic 
compliance to another orientation can fully populate the matrix. For 
example, rotation about the 3 axis gives shear/normal coupling, 
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Lekhnitski coefficients and variations in the rotated moduli E’, G' and 
v' aa ahown in Figure 2. 

The inelaatic reaponae of aingle cryatal materiala ia quite 
different from the behavior of polycryatalline nickel base auperalloys. 
The yield atreaa of aingle cryatal alloya ia a function of the material 
orientation relative to the direction of the applied atreaa. They also 
exhibit a aignificant tenaion/compreaaion aaymmetry in yielding. Figure 

3. [15] ahowa the yield atreaa of PWA 1 480 at 593°C aa a function of 
orientation along the [001] - [011] boundary of the standard 
aterographic triangle for tension and compression. The 
tenaion/compreaaion aaymmetry ia negligible near the [111] orientation. 
Further, the orientation dependence and tension/ compression asymmetry 
decrease aa a function of increasing temperature above a critical 
temperature. Aa an example, the variation of yield atreaa for PWA 1480 
at three orientations in tension and compression for a range of 
temperatures is ahown in Figure 4, [15]. In creep, at lower strain 
rates, aingle crystal alloya also exhibit orientation dependence even 
though the deformation mechanisms are expected to be different. 

Other observed behavior is similar to that seen in isotropic 
nickel baae auperalloys. These . materiala exhibit both strain rate 
sensitivity and cyclic hardening. Shown in Figure 5, [17] is the cyclic 

atreaa strain response for PWA 1480, 871°C, at different strain rates in 
the [111] orientation. An initial stress strain hysteresis loop as 
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FIGURE 2 - VARIATION OF ELASTIC MODULI WITH ORIENTATION FOR A .CUBIC SYMMETRIC MATERIAL 
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FIGURE 3. YIELD STRENGTH OF PWA 1480 at 593°C (1100°F) WITH 

ORIENTATION ALONG THE [001] - [Oil] BOUNDARY OF THE 
STANDARD STEREOGRAPHIC TRIANGLE. (REFERENCE [15]) 



TEMPERATURE, *C 



FIGURE 4. YIELD STRENGTH OF PWA 1480 AS A FUNCTION OF TEMPERATURE FOR 
THE [001], [Oil] and [111] ORIENTATIONS IN TENSION AND 
COMPRESSION. (REFERENCE [15]) 
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FIGURE 5. STRAIN RATE SENSITIVITY TEST OF PWA 1480 at 871 C FOR THE l 111] ORIENTATION. 

(REFERENCE [17]) 



compared to a loop after cyclic hardening is shown in Figure 6, [18], 

for PWA 1480 at 760°C in the [123] orientation. 

1 . 3 Constitutive Modeling History 

1.3.1 Continuum Mechanics Approach 

A number of investigators have proposed anisotropic plasticity and 
creep theories. The majority of these are based on a mathematical form 
similar to yield surface plasticity theory for isotropic materials, [19] 
■ [29]. Material anisotropy is introduced by including an anisotropy 
matrix in the yield function and kinematic hardening is modeled by 
utilizing a back stress tensor or a displacement of the yield surface. 
Isotropic hardening is accounted for with a change in yield surface 
size. A recent modification of this approach involves using a unified 
constitutive equation while retaining the basic mathematical structure 
for the flow law but without an explicit yield surface [38], [49], in 
order to demonstrate the general features of the approach it is useful 
to summarize one of these theories as an example. 

Lee, Zaverl, and Shih have developed an anisotropic plasticity and 
creep theory that successfully predicts many of the features of the 
inelastic response of Zircaloy, [23], [24], [26], and [30]. The yield 
function is an extension of Hill’s work [28], and is given by 

M i j ( a i —a i ) (C j-Oj) - k 2 - 0 (1.2) 


in a six dimensional stress space, a^. The anisotropy matrix 

describes the variation of yield stress with orientation. The parameter 
k defines an effective yield surface size and the back stress a. 
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FIGURE 6. CYCLIC HARDENING RESPONSE OF PWA 1480 IN THE [123] ORIENTATION AT 760°C. (REFERENCE 18) 


describes the initial tension/ compression strength differential and 
subsequent kinematic hardening. All of the parameters, k, c^, and M 


are functions of the plastic strain increments. The flow rule. 




(1.3) 


is a generalization of the Prandtl-Reuss equations for isotropic 
materials. The effective stress c is defined by 

• W'-i) ( v a .t ) 

and the effective plastic strain increment is defined by 

(d = P)2 ■ V E i dE j ( ’- 5) 

* 


where is the generalized inverse of and the generalized 

devlatoric stress is given by 


S i “ 3 M iJ ( °J V * (1 * 6) 

To extend the theory into the creep regime an effective creep 


— c — 

strain rate e is assumed to be a function of the effective stress o, 

defined by Equation (1.4), as well as the temperature and shear 

resistance t. The shear resistance is assumed to be proportional to k. 

The flow Equation (1.3) is modified slightly to become 


.C 

e i 



(1.7) 
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c 

where e are the creep strain rate components and the total inelastic 

strain is then the sum of the creep and plastic strain components. A 
method for the evaluation of material parameters has been established 
and the calculated and experimental response of Zircaloy 2 has been 
presented in References [23] and [26]. The determination of the yield 
function behavior requires load reversal tests and monotonic tension and 
compression tests in the three principal directions of anisotropy. The 
rate dependent material parameters are evaluated from load relaxation 
tests. 

A unified model that falls within the continuum mechanics approach 
was proposed by Stouffer and Bodner [38], for initial or deformation 
Induced anisotropy. The approach, as summarized in Appendix A, was 
investigated as part of the study for application to Rene N4. The model 
is based on transforming the anisotropic flow equation, when written in 
a six dimensional 3pace, into the eigenspace of the compliance (or 
stiffness) matrix. This results in a system of six uncoupled scalar 
flow equations in the compliance eigenvalues which has some advantage 
for computation. The major disadvantage is that the transf ormati on 
between the eigenspace and physical space depends upon the evolution of 
material properties arising from the deformation and is generally 
unknown. However, as presented in Appendix A, it is possible to 
determine this transformation exactly for various symmetry classes 
including cubic symmetry. The results in Appendix A show that once a 
specific transformation is established; for example, for cubic symmetry, 
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the response will be forced to follow this symmetry class exactly for 
all future deformations. Unfortunately, this condition appears to be 
too restrictive. Relaxation of the symmetry restriction leads to 
parameters in the eigenspace that are not physically motivated and 
difficult to determine experimentally. Thus, it was decided to abandon 
this eigenspace approach in favor of an approach with material 
parameters that are motivated by the physics of the deformation and 
directly related to experiments. 

The most significant shortcoming of the continuum mechanics 
approach as described in both of the above examples is that the actual 
deformation mechanisms are included in the theory only very crudely if 
at all. This limits the predictive capability of the model outside the 
range for which it has been calibrated. Another serious problem is 
modeling the plasticity, creep and fatigue interactions. Lee, Zaverl, 
and Shlh have made an attempt to correct this deficiency by assuming a 
relationship between the shear resistance in the creep equation and the 
yield surface size in the yield function. Finally, although not 
discussed above, the current theories do not contain methods of 
including temperature history effects. Current technology is based on 
evaluating model parameters from isothermal data at different 
temperatures and then interpolating between temperatures to determine 
the response. This does not account for a change in deformation 
mechanism due to a temperature excursion. 

The principal argument in favor of a continuum mechanics approach 
is the relative numerical simplicity. For the purpose of implementation 
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in a finite element code it is somewhat less complicated and generally 
requires less calculation than a crystallograpic approach. 

1.3.2 Crystallographic Approach 

Early developments in this approach are attributed to Taylor [31]. 
Bishop and Hill [32], [331 and Bishop [34]. The application of the 
crystallographic approach to single crystal nickel base superalloys 
began with the work of Paslay, Wells, Leverant and Burck [35], [36], and 
more recently by Shah [8] to the Y* phase of these alloys. 

The crystallographic approach is based on identifying the active 
slip planes and slip directions. The shear stresses are computed on 
each of the slip planes from the applied stress. The slip deformation 
is computed on each slip system and the macroscopic inelastic strain 
rates or strain Increments are then the sum of the contributions of the 
Individual slip systems. This approach is computationally intensive 
since there are a number of slip systems to be considered at each point 
in the body. Furthermore the response at the crystallographic level is 
not necessarily easy to determine. 

The classical constitutive assumption is Schmid's law in which it 
is assumed that the slip on a particular system is a function of the 
resolved shear stress on the slip plane in the direction of slip. This 
assumption, however, does not apply for LI ^ crystals in general or for 

single crystal nickel base 3uperalloys. More recently, Pope and others 
[7]-[12], [14]-[16], have proposed an extension to Schmid's law that 
includes three components of stress. This approach, which is discussed 
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in Chapter 2, appears to be successful in correlating the deformation 
mechanism with several of the macroscopic anisotropic effects. 

The principal advantage of this method is that a significant 
portion of the model is based on the. physics of the deformation 
mechanisms. Presumably, this will enhance the predictive capability of 
the model. Furthermore, as additional information is obtained about 
deformation mechanisms at different temperatures and strain rates the 
local constitutive models can be modified to accommodate the new 
knowledge. -The major disadvantage of this approach is that knowledge of 
the metallurgy and the Interface with mechanics must be understood. A 
less important objection is the additional difficulty in numerically 
implementing crys tallogr aphic models in finite element codes and the 
Increased computational requirements. 

1 . 4 Scope of the Present Study 

The purpose of this research is to provide a tool for the 
mechanical design and analysis of single crystal turbine blades and 
vanes. This objective is achieved through the development of an 
appropriate constitutive model and its implementation in a general 
purpose three dimensional nonlinear finite element code. 

The important features of the nonlinear constitutive response are 
to be modeled. These include the orientation dependent stress/strain 
behavior and tens ion/ compress ion asymmetry. The constitutive model must 
also predict cyclic and creep response. Because of the extremely 
limited data base, some important aspects of the material response have 
not been modeled in this study. The most important aspects of the 
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response that need further Investigation are the nonisothermal behavior 
and the response of the material in a multiaxial stress state. Also of 
concern and not currently included in the constitutive model is the 
cyclic softening/hardening behavior and latent hardening due to 
intersecting slip systems. 

The constitutive model has been developed using the 
crystallographic approach. The orientation and tension/ compression 
yield asymmetry can be accounted for using the metallurgical models 
developed in References [7]-[12], [14]-[16]. The development of a 
constitutive model from metallurgical concepts should provide a basis 
for further development of both the mechanics and metallurgy associated 
with single crystal technology. 

The intent in developing a finite element program is to produce a 
code that is general enough to model a wide range of geometries and load 
histories, flexible enough to easily modify the local constitutive 
equations and at the same time efficient. The need to model general 
three dimensional geometries lead to the choice of the twenty noded 
isoparametric solid element. The code efficiently and accurately 
integrates the constitutive equations over piecewise linear load 
histories using an initial strain method with a dynamic time 
Incrementing procedure. Because of the intended application, small 
displacement and small strain measures are utilized. 
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CHAPTER II 


DEVELOPMENT OF AN ANISOTROPIC INELASTIC CONSTITUTIVE EQUATION 
The philosophy in developing the constitutive model is to produce 
a system of equations that is as simple as possible while still 
describing the important aspects of the observed material response. As 
a result of an extremely limited data base, certain effects that could 
have been included are omitted rather than engaging in speculation. The 
effort in this report is primarily directed toward modeling Isothermal 
tensile and creep response with some hypothetical extensions to strain 
rate effects, stress relaxation and initial cyclic response. 

2. 1 Observed Deformation Characteristics 

Slip trace studies of single crystal alloys indicate that one or 
more types of slip may occur under different temperature, orientation and 
strain rate conditions. These Include: A. (Octahedral Slip) slip on the 

four octahedral planes in the three directions similar to the [10 1] 
direction (see Figure B2, Appendix B) ; B. (Octahedral Slip) 3lip on the 

octahedral planes in the three directions similar to the [1 2 f ] 
direction; and C. (Cube Slip) slip on the three cube planes in the two 

directions similar to the [T 0 1 ] direction. The slip conditions 
occurring during creep and tensile tests are examined for use in 
development of the model. 

2.1.1 Tensile Response 

Recent metallurgical research in the behavior of nickel base 
single crystal superalloys indicates that the tensile response of these 
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alloys is controlled to a large extent by the behavior of the Y ' phase 
alone, [7], [15] and [163- Thus, the response of L 12 ordered alloys and 

Y/Y' single crystals are discussed together. 

In uniaxial tests the critical resolved shear stress (CRSS) , or 
component of stress on the slip plane in the direction of slip required 
for yielding, is approximately constant or Increases slightly up to a 
critical temperature and is a function of orientation. Below this 

critical temperature slip occurs primarily on the octahedral (111)[1 0 
1] slip systems and significant orientation dependence and 
tension/compression asymmetry is observed. Above the critical 
temperature there is a sharp drop in the CRSS similar to the yield 

stress in Figure 4, cube slip, (010)[1 0 1] becomes more prominent, and 
the tens ion/ compress ion asymmetry is reduced, [16]. There are two 
exceptions to this behavior. First, test specimens near the [ill] 
orientation exhibit cube slip at all temperatures. Second, loading near 
the [0 0 1] orientation is thought to produce only octahedral slip at 
all temperatures since the resolved shear stress on the cube planes is 
zero. 

2.1.2 Creep Response 

The understanding of the creep properties and the active slip 
systems during creep in Y/Y' alloys is much less complete than the 
tensile response. Only a limited number of studies have been reported 
in the literature; however, a few important observations can be included 
in this review. 
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Unlike the tensile yield response, the creep properties of ordered 
L 12 crystals are much different than the creep response of the Y/Y ' 

alloys [693* The creep rate of Y* is significantly higher than that of 
Y/Y' alloys under identical conditions. 

Some nickel base superalloy single crystals exhibit an incubation 
period in addition to the usual primary, secondary and tertiary creep, 
see Figure 7b [7]. The incubation period and primary creep are 

attributed to (111) [T 2 7] slip. In single crystal Rene N4 at 

temperatures from 760°F to 1150°C the response did not include an 
incubation period similar to that in Mar-M200, see Figure 7 and 11. 

In secondary creep at temperatures above the critical temperature 

slip is inferred to occur by (111) [T 0 1] slip [8]. However, cube slip 
is considered to be present near the [111] orientation and absent near 
the [0 0 1] orientation similar to the tensile response. These 
deformation modes have been verified by transmission electron microscopy 
(TEM) studies for PWA 1 444, [63] and in addition dodecahedral slip, 
(011) planes, was observed for loading in the (001) and (Oil) 
directions. In this study a tens ion/ compress ion asymmetry was observed 
and typical response showed a creep rate that is higher in tension than 
compression for a creep in the [1 1 0] direction. This response is 
similar to the tension/corapr ession asymmetry observed in the yield 
stress. 
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FIGURE 7. CREEP CURVES OF MAR-M200 CRYSTALS AT 760°C AND 690 MPa. 

(REFERENCE 7) 


21 



2.2 Metallurgical Models 


A major step in understanding this behavior was made by Takeuchi 
and Kuramoto [10] in their study of Ni^Ga single crystals. They 

proposed that the increase in CRSS up to the critical temperature and 
non Schmid's law behavior is a result of cross slip of screw dislocation 
segments from the octahedral to cube planes. The cross slipped segments 
pin the dislocations and therefore increase the flow stress. They 
further proposed that the cross slip mechanism i3 thermally activated 

and is driven by the resolved shear stress in the (01 0) [ 1 0 1] system. 
However, their model did not fully, explain the tension/compression 
asymmetry or the orientation dependence observed in other L1 2 alloys. 

An improvement in the Takeuchi-Kuramoto model was proposed by 

Lall, Chin and Pope in [11]. In their theory the octahedral -| [T 0 1 ] 

dislocation is an extended dislocation consisting of two Shockley 

partial dislocation pairs, -g [ 2 1 1 ] + -g [ 1 1 2]. In order to cross 

slip the pair must constrict into a single ^[101] dislocation. The 

constriction is aided by a shear stress on the (111) plane in the 

[12 1] direction. It is important to note that a shear stress in the 
opposite direction extends the dislocation pair and thus inhibits cross 
slip. This effect is generally referred to as the "core width effect" 
and gives rise to the tension/ compression asymmetry observed in these 
alloys. Lall, Chin and Pope also proposed that the change in flow 
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stress on the octahedral plane in the [1 0 1] direction, for 

example, from a reference state at 0°K is given by 

4T ()n) -Aexp . (2.1) 

where A is a constant, T is temperature, and k is Boltzmann’s constant. 
The parameter H is a function of the resolved shear stress on the (ill) 

plane in the [121] direction, t 1 , and the resolved shear stress on the 

(010) plane in the [T 0 1 ] direction, x^. The stress components x 1 and 
are shown in Figures B1 and B2. By expanding H in a Taylor series 
about the reference condition H q , Equation (2.1) becomes 

_ H n + V i T i + V ? T ? + ••• 

At ( 111) - A exp 2 L-L , (2.2) 

where V 1 and are constants. 

In the orientations and at temperatures where cube slip is 
dominant the orientation dependence and tension/ compression asymmetry is 
reduced considerably. Therefore, it is believed that Schmid’s law can 
be used to relate the slip rate to the resolved shear stress on the cube 
planes. 

A major effort in this work is to recast the above results, both 
experimental and theoretical, into the form of a constitutive equation 
involving state variables that characterize the change in material 
properties due to the deformation history. These results are valuable 
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for determining the appropriate independent variables and functional 
forms . 

2. 3 Kinematics 

It is necessary to develop the kinematic equations that relate the 
shear stresses on the crys tallograpic planes and in the directions of 
interest to the stress tensor relative to the principal material axes. 
The constitutive equations are applied at the crystallographic level and 
it is therefore necessary to also develop the relationship between the 
crystallographic shear strain rates and the -macroscopic inelastic strain 
rate tensor. The required stress and infinitesimal strain relations 
were derived by Bishop in [3*0 and are summarized below. The sign 
convention used in this research and the details of the calculations are 
presented in Appendix B. 

The relative rate of displacement of a point located at x^ due to 

* a8 

a uniform shear strain rate Y parallel to the a plane in the 8 
direction is given by 


u 


i 


i 

J 


j J*i 


( 2 . 3 ) 


where n° is the unit normal to the a plane and is the unit vector in 

the 8 direction on the a plane (no summation of repeated indices). 
Summing for all planes and slip directions gives 


r ' a 8 a 0 a8 

u ‘ ' Jj T "j x j *i 


( 2 . 4 ) 


Applying the linear strain displacement relation 
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1 

£ 1J - 2 


(2.5) 
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gives the desired relationship between the crystallographic shear strain 


• gg 

rates Y , and global strain rates , 


. r 1 *aS,»a0 a .aB_cu 
e lJ " l 2 Y (t i n j + n i } 

ap 


( 2 . 6 ) 


The stress vector Oj on a plane a is given by 


°j • [ "u n l 


(2.7) 


where is the applied stress tensor in the principal material axes. 

The local shear 3tress component in the B direction on the a plane is 
given by 


- a8 I o.l a8 


T a 


J 


J J 


( 2 . 8 ) 


Combining the above two equations gives 


afi r _a .a 

T “ J a ij n i 


a „aB 


(2.9) 


the transformation between the global applied and local stress 
components. 

In summary, the constitutive model is based on separating the 
total global strain at a point in a body into elastic and inelastic 
components. The elastic strains are computed using Equation (1.1) and 
the applied stress tensor. Inelastic strains are obtained by 
integrating the inelastic strain rate with time. The inelastic strain 
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rate is calculated by summing the contributions from each of the 
crystallographic slip systems using Equation (2.6). The inelastic shear 
strain rate on each of the slip systems is computed from the inelastic 
constitutive equations which are developed in the remainder of this 
chapter. In addition, the inelastic shear strain rates on the active 
slip systems are a function of the local shear stresses which are 
computed from Equation (2.9). 

2.4 Inelastic Constitutive Equations 

There are at least two separate flow equations required to 
describe the tensile response of single crystals, i.e., cube slip and 

octrahedral slip in the [10 1] direction. In addition dodecahedral 

slip and octahedral slip in the [T T 2] direction have been reported in 
creep for some single crystal superalloys. Since no incubation period 
has been observed for Rene N4 in creep at any temperature, it appears 

that octahedral slip in the [T T 2] directions may not be present. Creep 

deformation is attributed to the immediate operation of (ill) [T 0 T] or 

(001 ) C 1 0 1] slip depending on orientation. Further, dodecahedral slip 
has not been observed for Rene N4 in a recent study of tensile and 
fatigue response, [6], [37]. Thus in this work it is assumed that only 

cube slip and octahedral slip in the [101] directions are active, and 
the model will be developed on this basis. 

In addition, there appear to be at least two deformation 
mechanisms present. At the higher strain rates, such as in tensile 
tests, deformation occurs by dislocations cutting the Y? precipitate. 
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The deformation rate depends on the size and volume fraction of Y'. The 
orientation dependence of the CRSS at temperatures below the critical 
temperature depends on the core width effect and the cross slip 
mechanism. The creep strength of Y/Y* alloys appears to be controlled 
by the Y/Y' interfaces in the material. This explains the difference in 
creep behavior of Y* and Y/Y’ alloys. The activation energy during 
creep is much higher than expected from a vacancy diffusion process. 
Tien, et al [70], proposed that during steady state creep the presence 
of Y* particles leads to a dislocation network that exhausts the vacancy 
concentration. Deformation occurs by the less favorable process of 
emitting and diffusing interstitials, which is consistent with the 
higher activation energies. This model has been supported by many 
observations as summarized in [7]. 

The modeling effort should include both the tensile and creep 
mechanisms. In addition, it appears there should be a strong 
interaction effect since the dislocation network and Y • particles are 
Important in both cases. Unfortunately, there is very little or no 
evidence available for the development of a model to characterize these 
interaction effects. In the next two subsections a constitutive 
equation for cube and octahedral slip is proposed. 

2.^.1 Octahedral Flow Equation 

The octahedral slip system is active for both tensile and creep 
deformation, but the deformation mechanisms are different as discussed 
above. The model for octahedral slip is proposed as the sum of two 
components i.e., 

27 



( 2 . 10 ) 



The first component is used to characterize dislocation cutting of the 
Y* particles and the second is motivated by the interstitial emission 
and diffusion mechanism. The two terms must be coupled to characterize 
the effect of the dislocation network on both mechanisms. Equation 
(2.10) is typical of classical models involving the use of both plastic 
and creep strain components. The origin of this approach is clearly 
based on the physical motivation of two or more deformation mechanisms. 
The thrust of the unified strain measure is to eliminate the need to 
identify the coupling between the two mechanisms. Unfortunately, many 
of the unified models are reasonable for modeling plasticity or creep 
but are not completely adequate for both. In the present study a 

* ciS 

coupling parameter is proposed in combination with the flow laws for y“ 
and ^2 that appears to be an improvement. 

The functional form of the octahedral flow equations for y“ ct is 

similar to the exponential form developed by Bodner et. al. for 
Isotropic materials. This form has been used with success for a number 
of materials including isotropic nickel base superalloys at high 
temperatures, [39], [40],' [50] = [60]. The Bodner flow equation was taken 
from an expression for mobile dislocation velocity proposed by Gilman, 
[61], [62] and is expected to capture the essential features of the 


28 



response for Rene N4 at the higher strain rates below the critical 
temperature. The form of the equation used in this study is 


- D. exp 
OCT . 1 K 


1 


zf \ ", 


a8 1 


a8 


aS i 


( 2 . 11 ) 


,a8 


where is a state variable Introduced to characterize the resistance 


to Inelastic flow. The parameters and n^ are used to characterize 


the limiting strain rate and the strain rate sensitivity. The 

octahedral flow equation and state variable must also characterize 

the orientation and tension/compression asymmetry properties. Since the 
form of the flow equation is similar to the Lall, Chin and Pope model 
their results can easily be incorporated into Equation (2.11). 


The functional form of the octahedral flow equation for y“ ct is 
taken in a form similar to Equation (2.11), l.e. 


■ D 2 


aB i 


aB 


|t“ b I 


( 2 . 12 ) 


This form is selected because it is consistent with representations used 
to characterize thermally driven diffusion processes as well as mobile 
dislocation velocities. Further, maintaining a similar structure allows 

for the same Interpretation of the state variables Z and Z°®,which 

should aid in establishing coupling between the two terms. 

Combining Equations (2.10), (2.11) and (2.12) gives 
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4 -1 -1 
The constants D 1 and D 2 are chosen to be 10 sec and 1 sec to allow 


uncoupling at high and low strain rates. The first term is negligible 
during creep and the second term is negligible during high rate tensile 
tests. Both terms are active at intermediate values of stress or strain 
rate. The constants n 1 and n 2 characterize the strain rate sensitivity 

and are determined independently. This is consistent with the observed 

response of single crystal alloys [7]. The state variables Z° 8 and Z° 8 

include work hardening arising from the development of a dislocation 
microstructure and include core width and cross slip effects. 

rc ft 

The flow resistance at high strain rates, Z 1 , is assumed to be 


similar to the Lall, Chin and Pope result, Equation (2.2), and is given 
by 


aft _aB ^ aB . ,, aB 

, - Z * ♦, - V,, t, * V 12 t 2 


( 2 . 1 ») 


ft rr ft 

where «j> 1 is the initial value of Z 1 , Z (0) - 0 and and V 12 are 

constants. The parameter Z a8 is a measure of work hardening and and 
aB 

t 2 are the shear stress components associated with the core width 
effect and the cross slip mechanism previously discussed. The shear 

ft (*Y ft 

stress t i is on the same octahedral plane and is perpendicular to t 
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The sign convention used for -r aB and t° B is given in Appendix B. The 

n ft 

shear stress is the magnitude of the resolved shear stress on the 

Q B 

cube plane in the same direction as -t . The increase in flow 

flft 

resistance due to work hardening Z , is given by the evolution equation 

Z aB - M 1 (W 1 -Z aB )T° B Y a8 ‘ (2.15) 

where M 1 and are material constants. The functional form of the 

ftfl * rtft 

evolution equation involving inelastic work rate t r rather than 

y alone is very similar to the form used in the Bodner model for 
isotropic materials. A major difference is that there is no recovery 
term Included in Equation (2.15). Also absent is an explicit function 
to model the effect of hardening due to intersecting slip systems. Even 
though some theoretical work has been done on latent hardening due to 
intersecting slip systems, [64]-[67], data to support the inclusion in 
the constitutive model currently do not exist for Rene N4. 

The expression for the state variable z“ at low strain rates is 

similar in form to Equation (2.14) and is given by 


,a8. 


aZ aB + $ 2 - 


V x aB V T a8 
21 T 1 + 22 t 2 


✓3J. 




( 2 , 16 ) 


31 



where 


and V 22 are material parameters associated with the shear 

stress components , and and J 2 is the second invariant of the 
deviatoric stress tensor. Equations (2.14) and (2.16) are coupled since 

a ft 

the same work hardening term Z appears in both equations. The 
parameters a and $ 2 define the development of the dislocation structure 

(hardening) at the stress-strain rate levels associated with creep. 
Equation (2.16) is valid for primary and secondary creep only. A 
thermal recovery term may also be necessary at high temperature. 
Tertiary creep is neglected because it occurs at strain levels beyond 
the range of interest for this study. 

It is possible to include a component of kinematic hardening into 

Cl 3 QL 3 

the model. Two values are assigned to each state variable, Z + and Z_ , 
corresponding to the direction of slip. Equation (2.15) i3 used to 


" s 

calculate z“ in the direction of slip. In the opposite direction the 


hardening rate, Z QD , is given by 



n 7 a& 

q Z SD 


(2.17) 


where q is a material parameter. Isotropic hardening occurs when q = 1 
and kinematic hardening occurs when q ■ - 1. This form is similar to 
the system used by Bodner et al [38] for tension/ compression response of 
isotropic materials. 
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2.4.2 Cube Flow Equation 


Cube slip occurs in specimens where the loading is near the 
[111] material direction and becomes increasingly important at high 
temperature. Furthermore, Schmid’s law is a good approximation for cube 
slip. 


• (jft 

The inelastic shear strain rate, on an a cube plane in the 

B direction is formulated with two terms similar to Equation (2.10) for 
octahedral slip and is given by 
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Once again the constants and are the limiting strain rates and are 
4 -1 -1 

chosen to be 10 sec and 1 sec to allow separation of the two terms. 
The constants and n^ reflect the strain rate sensitivity of the 

material in cube slip, and Z°^ and Z° 8 are the state variables that 
include work hardening. 

The evolution equations that account for the work hardening are 
similar to that for octahedral slip except that the orientation factors 
are excluded. These can be summarized as 


,oB 


+ Z 


aB 


(2.19) 


and 


_a8 . _aB 

Zjj - $ 4 + bZ 


( 2 . 20 ) 


a g 

where the evolution of Z is given by 
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rt A 

where Z p (0) - 0 and and are determined from the high rate data. 
The parameters ^ and b relate the hardening in the high rate and low 


rate response for cube slip. 



CHAPTER III 


EXPERIMENTAL DATA AND MATERIAL CONSTANTS 

Moat of the experimental data available for Rene N*< are from 
tensile and creep tests with the specimens oriented in the [001] 
material direction. Unfortunately, this is not particularly useful when 
developing a constitutive model for a material that exhibits orientation 
dependence and tens ion/ compress ion asymmetry. The two temperatures at 

which the best data base for Rene N4 is currently available are 760°C 

and 980°C. In this study the temperature is limited to 760°C because 
the orientation dependence and tension/compression asymmetry are more 
predominant and the repeatability of the data is better. 

At 760°the stress-strain response for several orientations in both 
tension and compression, and fatigue response was determined by Gabb, 

Gayda, Miner, and Voigt [6], [37]. Also available at 760 °C are tensile 
stress-strain and creep data at three orientations for an earlier 
version of Rene N4, designated as VF317, Reference [68]. The slightly 
different chemistry for VF317 may be the source of differences in the 
observed response. The tensile response reported by Gabb, Gayda, Minor 
and Voigt for specimens oriented in the [001] was about -30 percent 
weaker than VF317 in tension. $^nce the response characteristics are 
significantly different for the two data sets it was not possible to 
develop a single set of material constants for the constitutive model. 

The Gabb, Gayda, Miner and Voigt data were used to develop the 
constants and te3t the model for octahedral and cube slip at high strain 
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rates in tension, compression and fatigue. The VF317 data were used to 
develop the constants for octahedral slip at high strain rates in 
tension and low stress levels in creep. 

3.1 Rene N4 Composition A Response 

The tensile data reported by Gabb, Gayda, Miner and Voigt, 
Reference [6], consist of seven tests at three orientations. The 
specimens had a cylindrical gage section 19 mm long by 4.7 mm in 
diameter. The tests were run at a constant crosshead rate with an 

-4 -1 

Initial strain rate of about 2 x 10 sec . Data recorded were load 
and crosshead displacement. Plastic strain was estimated using the 
offset from the elastic loading line and the specimen gage length. The 
resulting stress strain curves are shown in Figure 8. 

It was confirmed by TEM analysis that slip in the [001] and 
[Oil] specimens was in the octahedral system while cube slip was 

observed with orientations near the [Til] orientation. Tensile axis 
rotation was observed in the specimens tested to failure. The 
variability of the yield stress within a crystal was small; however, the 
variability between crystals was as large as 23$. This difference which 
is attributed to the orientation and tension/compression asymmetry is 
typical of other single crystals near this temperature. However, at 

980°C the orientation and asymmetry properties are much less important 
and Schmid's law appears to be applicable, Reference [6]. Stresses 
beyond yield increased initially and flattened out for the [001] and 

[T 1 2] specimens. The [0 11] specimen exhibited decreasing loads 
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TENSILE RESPONE OF RENE N4 AT 760 C. (REFERENCE 6) 


sharp serrations to about 5? strain. The serrated yielding, which is 
observed in PWA 1480 and Mar M200, is attributed to the operation of a 
small number of slip planes, typically one, two or four. 

The fatigue response at 760°C reported in Part II of Reference 
C37], consisted of the monotonic yield points for specimens at six 
orientations in tension and compression and first cycle hysteresis loops 
for three orientations. The specimens had a cylindrical gage section 
15mm long by 5mm in diameter. The control waveform was sinusoidal with 
a frequency of .1 Hz. 

The tensile and compressive yield stresses at several orientations 
are shown in Table 2. The [001] specimen had the greatest initial 
yield in tension and the [ 011 ] specimen was strongest in compression. 
Both specimens displayed significant tension/ compression asymmetry while 

the [T 1 1], [023] and [2 3 6] specimens displayed very little or 
none . 

The initial hysteresis loops for [0 0 1], [0 1 1 ] and [T 1 1] 

orientations are shown in Figure 9. The [0 1 1 ], [T 4 5 ], [2 3 6 ] 
and [02 3] specimens had serrated flow characteristics. All specimens 
displayed slight hardening with continued cycling which increased with 
increasing strain range but was generally less than 10?. For all tests 
the response stabilized well before half life. 
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FIGURE 9. TYPICAL FIRST CYCLE HYSTERESIS LOOPS FOR RENE N4 AT 760 C 

(REFERENCE 37) 


3.2 Rene N4 VF317 Response 

The tensile data for Rene Nil VF317 at 760°C consist of three 
specimens at three orientations. The specimens had a Codep coated 
cylindrical gage section 18mm long by 3.3mm in diameter. The tensile 

tests were run at a constant strain rate of 8.333 x 10 sec to 
inelastic strains of .2$. At .2? yield the exten 3 ometers were removed 
and the specimens were pulled to failure at a constant head rate of 
.021 1 7 mm/ sec. 

The response to approximately 1? strain, as shown in Figure 10, 
displays considerable orientation dependence. The stress increased 
about 30? between the .02? yield stress and ultimate stress for all 
specimens. Although no TEM work was done to determine active slip 
systems it is reasonable to assume that cube slip is not active for 
these orientations. 

Tensile creep rupture tests were performed on Rene N4 VF317 at 

760° for the same three orientations. The 3.175 mm diameter specimens 
were Codep coated. A total of seven creep tests at three orientations 
were conducted and the results are shown in Figure 11. Two of the test 
at 620 MPa were terminated prior to failure. The reduction in area 
ranged from 17? to 48% . 
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FIGURE 11. CREEP RESPONSE OF RENE N4 VF 317 AT 760 C 




Note that the specimens in the [0 342 940] orientation were 
stronger in the tensile tests and in creep than expected when compared 
to the other data set or to PWA 1480 (Figure 3). The specimens at this 
orientation were nearly as strong as the specimens in the [001] 
orientation. The reason for this anomaly is not known. 

3.3 Derivation of the Material Constants 

The octahedral slip constants are derived from tests where cube 
slip is not present, for example, the [0 0 1] and [011] orientations. 
Conversely, the cube slip constants are derived from tests where 
octahedral slip is not present i.e., [Ill] orientation. In addition, 
the first term in the octahedral or cube flow equation is negligible in 
creep and the second term is negligible at high strain rates. When a 
choice of orientations is available for developing constants the data 
from the [001] orientation is used since it is the primary loading 
direction in turbine blade applications. 

Since the constitutive model is developed on the individual slip 
systems It is essential to know the relationship between the applied 
stress tensor and local shear stresses and the relationship between 
local slip rates and the global strain rate tensor. To simplify this 
task the kinematic equations in Appendix B were implemented in a 
computer code. For example, the shear stress components from an applied 
stress of 100 MPa in each of four directions is shown in Table 3. The 
local stress in every slip direction on each of the octahedral and cube 
planes Is shown. The notation is described in Appendix B. The stresses 

in the (111) [T 2 7 ] system are used only for determining t 1 in the 
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model, t 2 is determined from the(001)[1 0 1] slip system and t 13 

determined for either the ( I'M ) [ 1 0 1] or (001 )[7 0 1] slip system 
dependening on the active system. 

For the purpose of developing material constants it is also 
Important to choose test orientations where the local stress and strain 
rate on each of the slip planes can be determined. For example, a 
load In the [001] direction will produce equal stresses on eight of 

the (111) [T 0 1] slip systems and no stress on the other (111) [T 0 1] 
systems or the cube systems. It is expected that slip on all eight 
systems will be the same. By comparison, a load on a [0 342 940] 

specimen will have shear stresses in all of the (ill) [T 0 1] slip 
directions. Further, it is expected that these shear stress will 
produce slip on more than one system and the 3lip rates are not expected 
to be equal. Thus, the tests in the [0 0 1], [0 1 1] and [111] are 
preferred to evaluating constants due to the uniformity of stress and 
strain rate in the octahedral and cube systems. Table 3 also shows why 
only octahedral slip is found near the [001] orientation (no stress on 
the cube planes) and only cube slip is found near the [111] direction 
(dominant stress on the cube planes). 

3.3.1 Flow Equation Constants 

The octahedral and cube slip systems are assumed to have limiting 

4 -1 

strain rates, and D^, of 10 sec in high strain rate tests and the 
limiting strain rate in creep, D^, is chosen to be 1 sec 1 . Neither 
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data set contains creep tests near the Cl 11] orientation that would 
activate cube slip. As a result. It is not possible to evaluate any of 
the constants associated with cube slip in creep and is set to zero 

4 -1 

for all of the numerical exercises. The value 10 sec has been 
accepted as a constant for many materials except under extreme loading 

rates. The value D 2 - 1 sec 1 was chosen since creep rates are about 

-4 

10 less than tensile rates. 

The constants n 1 , and n^ are all evaluated using essentially 

the same method. To compute n 1 two or more tensile tests at different 

strain rates with the specimens in the [001] orientation are necessary. 
Two or more creep tests with specimens in the [001] orientation and at 
different stress levels are required to determine n 2 . The evaluation of 

n^ can be accomplished with two or more tensile tests at different 

strain rates with specimens oriented in the [111] direction. The 
tensile tests are loaded under constant strain rate conditions to 
failure and the saturated values of stress are determined. The 
secondary creep rates are used from the creep tests. 

At high strain rates where the contribution of the creep rate term 
is negligible. Equation (2.11) may be rewritten as 
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ex B 

where Z“ is the fully saturated value at that strain rate. The [001] 

* ctB ctB 

orientation has eight equally active octahedral systems and and t 

are easily calculated from the applied stress and the measured inelastic 

ci B 

is plotted against Zn t the slope 

n ft 

of the line will be - n^ . Once n^ is known the saturated value of z“ 

can be calculated from Equation (3-1) or the ordinate of the line. 

The constants n^ and Z°® are calculated in a similar manner from 

the secondary creep rates while neglecting the contribution from the 
high strain rate term. The values for n^ and n^ and the saturated 

values Z ° 6 and Z° S are obtained from tensile and creep tests using [ 111 ] 

oriented specimens. The parameters n^ , n 2> and n^ are independent of 

orientation and are also assumed constant for constant temperature. The 
ciB 06 

values z“ , and Z 2 depend on orientation. 

Since the required tensile tests were not available for either of 
the current data sets, the value for n^ was estimated from an analysis 

of PWA-1480 data at 1600°F and 1800°F. The value for n^ was assumed to 
be the same. Data were available for the calculation of n 2 for the 
Rene N4 VF317 chemistry. Obviously the strain rate sensitivity 


strain rate. When In 
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predicted by the current set of material constants is at best an 
estimate. 

- 3.3.2 Evaluation of Orientation Factors 

The orientation dependence and tension/compression asymmetry of 
the material in octahedral slip is characterized by the constants V 


and V at tensile strain rates and and V 22 in creep. The constants 

can be evaluated from any three tensile tests at constant strain rate 
(or any three creep tests) as long as the orientations and/or 
tension/compression sense is different and cube slip is not involved. 
Best results, however, are obtained- if the orientations are not close 
together and both tensile and compressive data are used. Tests close 
together tend to magnify the experimental variability. In addition, it 
is very helpful if the active slip systems are equally stressed as 
explained earlier. An optimum set of tests is probably [001] in tension 
and compression and [Oil] in either tension or compression. This choice 
of orientation is well away from the cube slip regime, there is a large 
variation in t 1 and x^, and the model is calibrated for the orientation 

of greatest interest. 

For the high strain rate tests where the contribution of the creep 
rate term is negligible the saturated value for for the three tests 


can be calculated. Rearranging Equation (2.11) gives 
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where and t 01 ^ are the values when the inelastic strain rate and 


stress become constant. Letting 


n ft 

S 1 - Z + $ - constant 


(3.3) 


Equation (2.14) can be written for each of the three tests, i.e. 


K 6 ), - - M'f), * v , 4 t “ s ), 
[<) 2 - s > *». 2 K 6 ) ; 


(3.4) 


KV s . ■ 

Since and are known from Table 3 or equivalent and z“ can be 

computed from Equation (3-2) for each test, the constants S , V 1 1 and 
V 12 can be determined. 

The calculation of the orientation factors and for the 

creep regime is similar. Neglecting the contribution of the high strain 
rate term. Equation (2.12) can be rearranged to give 
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(3.5) 


where y qqj is the value during secondary creep. Once again letting 


S 2 - aZ°^ + 4> 2 


(3.6) 


Equation (2.16) can be written for each of the three tests, i.e. 


50 




Since t?® and J 2 are known and can be computed from Equation 

(3.5) for each teat, the constants S 2 « V 21 , and V 22 can be calculated. 

For the Rene VF317 data set, the tensile stress strain curves 
for the three orientations were used to obtain the constants V 1 1 , 

and ^ 12 . These tests were not optimal for this calculation since the 

strain rate was changed before the stress level was saturated and none 
of the tests were in compression. The constants S 2# V 21 , and V 22 were 

calculated for the secondary creep rates for three specimens in 
different orientations. These tests also were not optimal since none 
were in compression. 

The Rene N4 data of Gabb, Gayda, Miner and Voigt contains only two 
stress-strain curves at a constant rate that activated octahedral slip. 

Recall that the test in the [Tl2] orientation activated cube slip. This 
is not sufficient to establish the orientation factors. Alternatively, 
the data set contains yield stresses for the initial quarter cycle of 
the fatigue loops for specimens in several orientations. Since the 
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control wave form is sinusoidal, the total strain range varies from test 
to test and the strain rate varies during the loading quarter cycle, it 
is not possible to determine the orientation factors as described above. 

Thus, to test the orientation and tension/ compression aspects of 
the constitutive model, some assumptions were made. It was assumed that 
the yield stress variation in the sinusoidal first quarter cycle is 
representative of the saturated stress in constant strain rate tensile 

-4 

tests. Furthermore, a strain rate of 2 x 10 was assumed. The 
constants were calculated from [011] tension and [001] tension and 
compression yield stresses from the first cycle. It is obvious that 
this procedure will not generate the actual constants, but it does 
provide a means for exercising the constitutive model with the data that 
are available. 

3.3.3 Evaluation of the Hardening Parameters 

The evaluation of the constants in the state variable evolution 
Equations, (2.15) and (2.21), for both the octahedral and cube slip 
systems proceeds in a similar way. The evaluation of M 1 , and (j> 1 for 

octahedral slip and and for cube slip i? based using the high 

strain rate tensile data. 

For the octahedral system. Equation (2.15) can be integrated to 

give 

Z aB - W 1 [1 - exp .(- Ml T aB Y aB )] . (3.8) 

Combining Equations (2.14), ( 3 . 3 ) and (3.8) yields 
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ln(W 1 ) - Ml T a6 T aB - in [S 1 - V 1lT ® 8 + V 12 t“ 8 - z“ 8 ] . (3-9) 

The quantity in [S 1 - V^t^ 8 ♦ V^t® 8 ~ is log l inear in the local 

rate of inelastic working, t Y . The parameter Z“ can be evaluated 

using Equation (3.2) and stress-strain test data since x , Y , D 1 and 
n 1 are all known. By plotting Equation (3.9) for the [001] tensile 
data. Figure 12, the constants and M 1 are obtained. Notice that at 

T a8 Y aB - 0, in (Wj) - ln[S 1 - V n x® 8 + V 12 t“ 8 " Z ? 8 J- 0nce W i 13 

obtained, <j> 1 can be calculated from Equation (3*3). 

In cube slip the calculation of M^, and Is identical to the 

procedure for octahedral slip. Since the creep strain rate is 
negligible at high strain rates the cube flow Equation (2.18), becomes 


*cx8 

Y - exp 




a8 1 


( 3 . 10 ) 


for a high 3 train rate test with the specimen loaded in the [ 111 ] 
orientation. Rearranging Equation (3.10) gives 


,b8 


- in 


tft\ 

l D 3 I 


3 x 


aS i 


(3.1D 


and at saturation 

,a8 


Z~“ - ^3 * W 3 " constant 


(3.12) 
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FIGURE 1 


Integrating the evolution Equation (2.21) results in 

Z® 8 - W 3 [1 - exp (- M 3 T aS Y a8 )] . (3.13) 

Using Equations (2.19) and (3-13) gives 

in(W 3 ) - M 3 -r a6 7 aB - In [Z® 8 - $ 3 ] . (3.1 *0 

r* A 

The parameter can be evaluated using Equation (3.11) and the stress 

ci 0 * ci B 

strain test data since t , t , D 3 and n 3 are all known. A plot 
similar to Figure 12 is used to find W 3 and M 3> Using the saturated 

n A 

value of z“ , Equation (3-12) is used to find $ 3 . 

Unfortunately, a constant strain rate tensile test with loading in 
the [111] direction was not available and M 3 and W 3 were assumed to be 

the same as for octahedral slip for the numerical exercises. 

3.3.4 Tensile/Creep Coupling Terms 

The tensile/ creep coupling constants, a and ^ for octahedral slip 

and b and ^ for cube slip can be calculated from the tensile constants 

and creep response. The constants a and $ 2 are calculated for an 

orientation where octahedral slip is dominant, such as [001]. The 
calculation of b and requires data for an orientation where cube slip 

is dominant, such as [111]. 

The value of z“ can be obtained from Equation (3.5) for the 
primary and secondary creep rates. At the beginning of the creep test 
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n ft 

prior to any hardening, the hardening parameter Z will be zero and 
Equation (2.16) reduces to 



V V 

21 t 1 + 22 t 2 


/3Jl 




(3.15) 


aS 


so <j> 2 can be calculated. During secondary creep Z is saturated and 


equal to W 1 , Equation (2.16) becomes 


if - aw, . * 2 


V T°^ V T a ® 
21 1 + 22 t 2 




viz. 


(3.16) 


and a may be calculated. 

The constants for cube slip in creep may be similarly calculated 
using the initial creep strain rate and the secondary creep strain rate. 
The data were not available in the current data base for the calculation 
of b and <fiy. 


A summary of all the constants for both data sets is presented in 
Table 4. There are a total of 22 constants in addition to the 3 elastic 
constants. However, D 1 , Dg. and are assumed known, thus there are 

effectively 18 constants required to model the inelastic response. 
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TABLE 4 - MATERIAL CONSTANTS 


CONSTANT 


GABB, GATDA, MINER 

VF317 VOIGT DATA 



10000 SEC 

2.02 

1 SEC" 1 

.708 

959.1 MPA 
3-7454 
3.5799 

5.656 MPA' 

177.9 Mpa 

52.52 
5626 MPa 

20070 MPa 

18690 MPa 


10000 SEC 

2.02 


1109 MPA 
.5752 
“.5799 

8.194 MPA' 
437.5 MPA 


10000 sec 

2.02 


1198 MPA 

8.194 MPa 
437.5 MPa 


The elastic constants at 760°C are E - 100.0 GPa, G - 96.5 GPa and 
v * .38. 
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CHAPTER IV 


FINITE ELEMENT IMPLEMENTATION 

The finite element implementation of the constitutive model was 
done using a twenty noded isoparametric solid element. This element was 
chosen because it is possible to model almost any three dimensional 
geometry and also allows for any orientation of the material principal 
axes. Order two Gaussian integration was used both for stiffness 
generation and the calculation of body forces. 

The ability to model piecewise linear load histories was also 
included in the finite element code. This capability is particularly 
useful when modeling stress strain tests or fatigue loops and also for 
certain analysis applications. Since the inelastic strain rate could be 
expected to change dramatically during a linear load history it is 
Important to include a dynamic time Incrementing procedure. The current 
study is restricted to isothermal conditions, consequently the finite 
element code is limited to steady- state thermal conditions during a load 
case. The term "load case" is used to denote a time period for which 
the initial and final loads and boundary conditions are defined and vary 
linearly between the end points. 

4. 1 Initial Strain Method 

The initial strain method is an efficient technique for 
incorporating time dependent constitutive models in nonlinear finite 
element codes [71], [72]. It is an economical technique because it is 
necessary to reform the stiffness matrix only for temperature changes. 
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The derivation of the nodal equilibrium equations follows directly 
from the principle of virtual work as described below. Element strains 
T 

{e }, in terms of nodal displacement {x} are given by 

{ e T } - CB]{x} , (4.1) 

where [B] is the element strain displacement matrix. Since the total 

£ 

strains are the sum of the elastic strains {e } and the inelastic 
strains {e 1 } it follows that 

{ e T ) - {e E } ♦ U 1 } . (4.2) 

The stresses are related to the elastic strains by 

la) - [E] {e E } . (4.3) 

where [£] is the matrix of elastic constants. Defining { f } as the 

A 

element nodal forces and (x} as the associated virtual displacements and 
applying the method of virtual work yields 

Ix) T tf} - / v {e} T {a} dV , (4.4) 

A A 

where {e} are the strains associated with the displacements { x } . The 
integration is over the volume of the element. Substituting Equations 
(4.1), (4.2) and (4.3) into (4.4) gives 

{x} T {f} - / v ([B]{x}) T [E]({e} - {e 1 }) dV . (4.5) 

A 

Since { x } is arbitrary and independent of the integration (4.5) becomes 
Cf} ♦ {f 1 } - Ck]{x} , (4.6) 

where the inelastic pseudo force {f^} is defined as 



and the elastic stiffness matrix [k] is defined by 

[Kr] - / y [B] T [E][B] dV . (4.8) 

When the elemental equations are assembled the global equilibrium 
equation 

[K]{d T } - {F} * {F 1 } (4.9) 

T 

is obtained. The matrix [K] is the global stiffness matrix, {d } 13 the 
total displacement vector, { F} is the vector of applied thermomechanical 

forces, and {F*} is the vector of inelastic pseudo forces. 

The calculation of the inelastic strain rates and the state 
variable evolution rate is accomplished in the constitutive subroutines. 
The integration over the element volume is generally performed 
numerically and in this study order two Gaussian quadrature was used. 
The total Inelastic force vector is then assembled by summing the 
contributions from all of the elements. 

4.2 Linear Load History 

In order to incorporate linear load histories into this scheme the 
total displacement vector is decomposed into elastic and inelastic 
components i.e., let 

Cd T } - ld E } ♦ {d 1 } . (4.10) 

g 

The vector {d ) is the displacement due to applied thermomechanical 
forces and {d*} is the displacement due to the inelastic pseudo-forces. 
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These displacement vectors can be calculated using 

(d E } - [K]" T {F} . (4.11) 

and 

{d*} - [K]” 1 {F 1 } . (4.12) 

The elastic displacements are obtained for the initial and final 
thermomechanical loads in the load case. Assuming a linear variation 
the elastic displacements at any time in the load case are given by 

{ d E } - {d E } ♦ (ll2_)({ d E } - { d E } ) . (4.13) 

O i O 

r o 

E '* E 

The vectors {d } and {d }_ are the elastic displacements due to the 

O r 

initial and final applied thermomechanical forces. The current time in 

the load case is t, and t and t. are the initial and final times in the 

o f 

load case. The displacements due to the inelastic strains at any time 
during the load case are given by 

{d 1 } - {d*} ♦ l {Ad 1 } . (4.14) 

TIME 

INCREMENTS 

The vector {d*} is the vector of displacements due to inelastic strains 
o 

at the beginning of the load case, and {Ad*} is a displacement increment 
due to the inelastic strains during a time step. The increment in 

displacements {Ad*} due to the change in inelastic strains {Ae*} during 
a time step is computed using 

{Ad*} - [K]” 1 { AF* } . (4.15) 
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The inelastic pseudo force increment {AF*} is calculated from 

UF 1 } - l / CB] T [E]{Ae : } dV , (4.16) 

ELEMENTS 

where {Ae*} is the change in inelastic strain during the time increment. 
4.3 Integration of the Constitutive Equations 

Although the integration of the constitutive equations to obtain 
the inelastic strain increment could be achieved by any number of 
schemes, [4.31. [4.4], the second order Adams-Moulton predictor- 
corrector method lends itself readily to the required equilibrium 
iteration. The global inelastic strain and the state variables (written 
as 2^ for convenience) are integrated using 


I ^i”^i-i ^ X I 

{Ae X } m - - f ( {e } j_.| ♦ Ce i ) i ) 


and 


(4.17) 


^ t i” t i-1^ 

Z . — ± — i-j— (Z + z ) + z 
i 2 ' 1-1 i' i-1 


(4.18) 


The vectors 


{e I } 1 and U* 1 }^ 


are the Inelastic strain rates at times 


t^ and t A-1 and Z ^ and Z 1 _ 1 are the state variable at times ^ and t 1 _ 1 . 

Since the Inelastic response of the material is computed on the 
crystallographic planes, there is a state variable for each slip system 
at every integration point. The inelastic strain rate is computed by 
summing the contribution from each of the active slip systems. The 
inelastic strain rate on each of the active slip systems is a function 
of the stress state and the state variable for each system as described 
in Chapter II. 
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4.4 Iteration Procedure 


At the beginning of a load case the initial and final elastic 
displacements are computed using Equation (4.11) and the displacements 
due to prior inelastic strains are computed using Equation (4.12). The 
total strains at the beginning of the load case are recovered for each 
integration point and the elastic strains are computed from 

{e E } - U T } - {e aAT } - {e 1 } , (4.19) 

where {e 1 } are the total strains and le° } are the thermal strains. 

The stresses are computed, transformed to the material axis system and 
Equation (2.9) is used to find the shear stresses on each of the slip 
systems. 

Using the current values of the state variables and the stress 
state allows the calculation of the Initial values of the shear strain 
rates and state variable evolution rates on each of the slip systems. 
Equation (2.6) is then used to compute the macroscopic inelastic strain 
rates. Before entering the time loop an initial time increment is 
computed and the inelastic strain increments are estimated using a 
forward Euler integration formula. From the estimated inelastic strain 
increments an Initial estimate is made for the inelastic pseudo force 
increment using Equation (4.16). The usual technique employed with the 
initial strain method is to assume that the incremental inelastic force 

{AF*}, the corresponding displacements {Ad*}, and che inelastic strain 

increments {Ae*} are all zero on the first iteration of a time step. 
The stability of the method can be improved considerably when a forward 
Euler integration of the Inelastic strain rates is used to make an 
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estimate of {Ac*}, {AF*}, and {Ad*} on the first iteration. This method 
results in an initial estimate which is much closer to the solution. In 
sample cases the overall number of iterations was reduced by more than 
one half. 

The procedure during a time increment is to estimate the solution 
on the first iteration using a forward Euler scheme as outlined above. 
Then displacements, strains, stresses, inelastic strain rates, and state 
variable evolution rates are computed at the end of the time increment. 
The inelastic strains and state variables are Integrated over the time 
increment and an improved inelastic force is computed. The procedure is 
repeated until convergence is achieved at the end of the time increment. 
Figure 13 summarizes the logic. 

4.5 Dynamic Time Incrementing 

In a computer code that allows a linear variation of loads with 
time a dynamic time Incrementing scheme is very desirable since large 
excursions in stress and inelastic strain rate are to be expected. The 
procedure used to compute the time increments requires a certain amount 
of initial experimentation to determine appropriate time step control 
parameters. However, once this has been done the procedure works quite 
well and is a tremendous improvement in economy over a constant time 
increment. 

The time increment is based on three separate time step control 
criteria. These are the maximum stress increment, maximum inelastic 
strain increment, and maximum rate of change of the Inelastic strain 
rate. The minimum time step calculated from the three criteria is the 
value used. Since the calculations are based on values readily 
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FIGURE 13. FLOW CHART OF FINITE ELEMENT SOLUTION PROCEDURE 
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available from the previous time step, little computational effort is 
required. 


4.5.1 Stress Increment Criterion 


A maximum stress increment criterion is used to control the time 
increment during primarily elastic excursions. This criterion is 
necessary to prevent overshoot of the point where significant inelastic 
strain rates begin. The calculation for the time increment is given by 

Ao, 




INC 
( A °k-1 } KAX 


(4.20) 


where At k-1 is the previous time increment, ( Ac k _ 1 ) MAX * 3 the maximum 


change in effective stress for all integration points during the 
previous time increment, and Ac^ NC is the maximum desired stress 

increment. The value for Ao INC is program input and will vary somewhat, 

depending on material constants. Typical values are about 15 MPa. 

4.5.2 Inelastic Strain Increment Criterion 

The inelastic strain increment criterion controls the time step 
when the stress and Inelastic strain rates are not changing 
significantly. This is given by 


At k 


Ac 


At, 


INC 


k-1 


(Ae k-1 } MAX 


(4.21) 


The maximum change in effective inelastic strain for 
points during the previous time increment is (Ae k _ 1 ) MAX 


all integration 
and Ae* _ is the 

INC 


maximum desired inelastic strain increment. The value for Ae? M „ is 

INC 

program input and typical values are about .000100. 
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4.5.3 Rate of Change of the Inelastic Strain Rate 


This criterion controls the time increment when the inelastic 
strain rate is changing rapidly such as in the "knee" of a 3tress strain 

curve. The quantity e 1 is a measure of how close the initial forward 
Euler estimation is to the final converged solution. The backward 
difference formula 


E i “ (e i-1 " E i-2 )/At k-1 


(4.22) 


"I “I 

is used to estimate e^ The maximum value of e for all Integration 


points (e )^^ 1,3 1136(1 t0 estimate the next time step using 



2At 


k-1 


Ac 


INC 


, ~I. 

' E i MAX 


(4.23) 


The parameter e is the maximum desired error by which the initial 
forward Euler estimation is in error. The value for e is program input 
and typical values are about .01. Equation (4.23) is derived simply 
from taking the difference between an Euler integration scheme and the 
more accurate second order Adams-Moulton method. 

4.6 Convergence Criteria 

Convergence is required at every Integration point. Two separate 
criteria must be satisfied for convergence. First, the difference in 
the effective inelastic strain increment from subsequent iterations must 
be less than a prescribed value, i.e. 

|AeJ - Aej^l < 6 e . (4.24) 
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Second the change in effective stress from subsequent iterations must 
also be less than a prescribed value, i.e. 



( 4 . 25 ) 


For all of the calculations done in this study the convergence 
tolerances were set at one percent of the maximum desired effective 
stress increment and one percent of the maximum desired inelastic strain 
increment used in computing the time step. 
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CHAPTER V 


COMPARISON OF THE CALCULATED AND EXPERIMENTAL RESULTS 
All of the numerical calculations were performed using the finite 
element code described in the previous chapter. Since all of the 
calculations were for uniaxial tests it was not necessary to model the 
actual specimen geometry. However, to test the force and stiffness 
assembly routines a model with two twenty noded bricks was used to 
simulate the specimen response. The calculations were performed to 
compare with test results and also to exercise the constitutive model 
and computer code under conditions for which no data are available. 

5.1 Comparison of Experimental and Calculated Response for Rene N 4 , 
Composition A 

The calculations were performed using the constants derived from 
the data presented by Gabb, Gayda, Miner and Voigt in References [6] and 
[37]. The constants are for octahedral and cube slip at high strain 
rate (Table 4). The features of the model that are exercised are the 
orientation dependence, tension/compression asymmetry, strain rate 
dependence and isotropic vs kinematic hardening for hysteresis loops. 
Because of the assumptions required in the derivation of the material 
constants due to the limited data base, most of the exercises for this 
data set are intended to demonstrate features of the constitutive model 
rather than correlation with experimental data. 
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5.1.1 Orientation Dependence and Tens Ion/ Compress Ion Asymmetry 

Recall that in Chapter III the constants were derived by 

hypothesizing that the .02 % yield stress for the first quarter cycle of 
sinusoidal loading was representative of the orientation dependence and 
tension/ compress ion asymmetry for constant strain rate tests. The 
octahedral constants were derived using the [001] tension and 
compression tests and the [011] tension data. The cube slip constants 

were derived using the [111] tension data. The stress response was 
calculated for these tests and for the tests at other orientations. The 
calculated results are compared to the experimental results in Table 5. 
All of the predicted results are within 6J of the experimental values 
and within the confidence level for the experiments. It can be 
concluded that the constitutive model successfully predicts the 
orientation dependence and tension/ compression asymmetry in octahedral 
slip. It can also be seen that the octahedral and cube flow equations 
produce good results when simultaneously activated. 

5.1.2 Cyclic Response 

The first 1.25 cycles of the fatigue response presented in 
Reference [37] are compared to the predicted results using both 
isotropic hardening and kinematic hardening (Figures 14 and 15). The 
stress ranges for the kinematic hardening assumption are generally 
better than for the isotropic hardening assumption. 
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TABLE 5. COMPARISON OF MONTONIC HELD STRESSES (REFERENCE [37]) 
AND CALCULATED SATURATION STRESSES WITH OCTAHEDRAL AND 

CUBE SLIP CONSTITUTIVE MODELS ACTIVE, RENE N4, 760°C 


TENSION OR .02* YIELD CALCULATED 

ORIENTATION COMPRESSION STRESS (REF [37]) SATURATION STRESS j ERROR | 
(MPa) (MPa) (X) 


[001] 

T 

956 

956 

0 

[001 ]* 

C 

>818 

-819 

.1 

[oil]* 

T 

748 

752 

.5 

[Oil] 

C 

>905 

-865 

4.4 

[Tii]* 

T 

817 

827 

1.3 

[Tn] 

C 

-842 

-828 

1.7 

[023] 

T 

695 

705 

1.5 

[023] 

C 

-747 

-741 

.9 

[236] 

T 

716 

725 

1.2 

[236] 

C 

-714 

-752 

5.3 

[145] 

T 

656 

692 

5.5 

[t«5] 

C 

-792 

-763 

3.6 

CONSTANTS WERE 

DERIVED 

USING THIS DATA 
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1 1-1000 1 

FIGURE 1A. PREDICTED AND EXPERIMENTAL RESPONSE OF FIRST CYCLE FATIGUE LOOPS FOR RENE NA AT 760°C 

ISOTROPIC HARDENING 


1000 t STRESS, MPa 
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KINEMATIC HARDENING 



There is a significant difference in the shapes of the curves. 
The use of a sinusoidal wave form has a tendency to produce flat curves 
but the predicted response even shows a stress drop near the peak strain 
points. Probably the most significant variation between the predicted 
and experimental results is seen in the response of the specimen in the 

[Til] orientation after the first 1/4 cycle. These differences can at 
least partially be attributed to the assumptions that were required to 
derive the constants. Recall, that the values for n 1 and n 2 were 

estimated from PWA 1480 data at 871°e and 982°c. This may not be 

representative of Rene N4 at 760°c. In addition, the hardening 
constants were estimated from Rene N4 VF317. These extrapolations could 
cause significant errors In the predicted hardening and strain rate 
sensitivity. 

A alight mismatch In elastic modulus is seen between the predicted 
and experimental results for the specimen oriented in the [011] 
direction. The elastic moduli used in the analysis are the nominal 
values and are entered in the principal material directions. A slight 
variation In elastic modulus from the nominal value is not surprising. 
5.1.3 Strain Rate Sensitivity 

In order to demonstrate the strain rate sensitivity of the 
constitutive model, first cycle fatigue loops were calculated for two 
different orientations and three different rates using the sinusoidal 
total strain history. The results of this exercise are shown in Figure 
16. Qualitatively, it can be stated that the observed trend is as 
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expected. Since no data are available, quantitative comparisons are not 
possible. 

5.2 Comparison of Experimental and Calculated Response for Rene N4 
VF317, Composition B 

The calculations were performed using the constants derived from 
the Rene N4 VF317 data. Reference [68]. The constants are for 
octahedral slip at high strain rate and creep (Table 4). The features 
of the model that are exercised are orientation dependent stress strain 
curves, orientation dependent creep response and stress relaxation 
behavior. 

5.2.1 Orientation Dependent Stress Strain Curves 

The stress strain response was calculated for a constant strain 

rate of 8.333 x 10 sec for three orientations. The calculated 
response compared to the experimental data is shown in Figure 17. 
The orientation constants were determined using all three curves so the 
variation in the yield stress levels is reproduced well. The hardening 
parameters were calculated from the [001] curve only. 

Notice that there is a significant mismatch in the elastic modulus 
for the predicted and experimental response for the [0 342 940] 
orientation. Recall the observed VF317 response does not agree with the 
PWA 1 480 response shown in Figure 3 because the yield stress for the 

[0 342 940] specimen, 20° from [001] in the [011] direction, is far 
above the expected level. The difference in the experimental and 
calculated elastic modulus suggests that the specimens in this data set 
designated as being in the [0 342 940] direction are probably much 
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FIGURE 17. PREDICTED AND EXPERIMENTAL STRESS STRAIN CURVES FOR 
RENE N4 VF317 AT 760°C. 
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closer to the [001] direction. This inconsistency is further verified 
through an example compression calculation. The predicted saturation 
stress response for a [001 ] specimen in compression was less than a 
third of the corresponding tensile [001] response. This is not 
realistic, especially since the model worked very well for the previous 
tension/ compression exercises 

5.2.2 Creep Response 

The tensile creep curves for seven specimens at three orientations 
were calculated and are compared with experimental data in Figures 18- 
20. The experimental data are from creep rupture tests which usually do 
not yield reliable data at the beginning of a test. Usually this type 
of test produces significant variation from test to test due to poor 
alignment and strain measuring systems as well as inherent differences 
from specimen to specimen. The correlation between the predicted and 
actual response for this type of test must be considered to be 
excellent. 

5.2.3 Stress Relaxation 

Although no stress relaxation data are available the stress 
relaxation response for specimens oriented in the [001] and [011] 
directions was calculated. Both predictions are for a constant 
displacement boundary condition with an initial stress level of 758 MPa. 
The qualitative behavior is as expected with the stress in the [011] 
oriented specimen relaxing faster than in the [001] oriented specimen, 
Figure 21 . The strain rates in the [001 ] oriented specimen are in the 
3ame range as in the creep tests. However, the strain rates in the 
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FIGURE 18. CREEP RESPONSE OF RENE M VF317 LOADED IN THE J001J DIRECTION AT 760 C. 
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FIGURE 19. CREEP RESPONSE OF RENE N4 VF317 LOADED IN THE [110] DIRECTION AT 760°C. 


FIGURE 20. CREEP RESPONSE OF RENE N4 VF317 LOADED IN THE [0 342 940] DIRECTION 

AT 760°C. 
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[Oil] oriented specimen span the range between the tensile and creep 
data. Although the model behaved as expected it would be interesting to 
compare it to actual data in this intermediate strain rate range. 
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CHAPTER VI 


DISCUSSION AND SUMMARY 

The constitutive model as developed in this report is physically 
based and is capable of modeling most of the relevant material behavior 
at constant temperature. However, consideration must be given to future 
development needs. In addition, it is useful to specify the data base 
required to fully characterize the material parameters in the 
constitutive model as it now exists. 

6, 1 Constitutive Model Development 

In the short term it should be possible to utilize the current 
constitutive model and finite element code as a viable mechanical 
analysis tool with little further development. For the future a number 
of potential refinements in the constitutive model are anticipated as 
outlined below. 

6.1.1 Coupling of Creep and High Strain Rate Response 
The work hardening in the flow equation for each slip system is 

__ fl 

characterized by a state variable Z . This state variable is a measure 
of flow resistance at high strain rates and in creep. In the current 
data base this approach appeared satisfactory, but questions remain for 
load histories in which both high strain rates and creep occur. To 
evaluate the coupling as it is formulated in this study requires a high 
strain rate loading followed by creep and/or a creep test followed by a 
tensile test. 

Since different deformation mechanisms are postulated for high 
strain rates and creep it is reasonable to assume that two work 


84 





hardening state variables may be necessary. If this is the case, it is 
logical that coupling would occur in the state variable evolution 
equations. 

6.1.2 Nonlsothermal Response 

The most obvious problem in modeling the noni sothermal response of 
Rene N4 is the lack of a data base sufficient to derive the material 
parameters over a range of temperatures. A simple nonlsothermal 
implementation can be accomplished in the same manner as an isotropic 
formulation. The material constants are interpolated for temperature, 
and temperature rate and history effects are neglected. 

Although the model has not been calibrated for any other 
temperatures it is possible to propose a method for interpolating some 
of the material parameters. In the flow Equations (2.13) and (2.18) the 
parameters D 1 , D^, and are Interpreted as limiting strain rates 

and are expected to be constant with temperature. The strain rate 
sensitivity parameters n 1 , n 2> n^ and n^ can be represented by an 

equation of the form 

n « | + b (6.1) 

as proposed by Bodner in Reference [53], where a and b are constants and 
T is the absolute temperature. The thermal dependence of the 
orientation parameters , <P 2> 1 , V 12 , V 21 and V 22 should be in the 

form of an Arhennius function, (Equation (2.2)). The variation of the 
material parameters in the state variable evolution Equations (2.15) and 
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(2.21) with temperature is unknown and would be postulated from the 
data. 

The validity of this approach needs to be investigated. A 
reasonable approach would be to compare predictions of the current model 
with nonisothermal experiments. It is anticipated that the state 
variable evolution Equations (2.15) and (2.21) may require a thermal 
rate term to Include thermal history effects. 

6.1.3 Thermal Recovery 

The treatment of secondary creep as a balance between strain 
hardening and recovery was originally proposed by Orowan in 1947, 
Reference [76]. Nearly all of the recently developed unified 
constitutive models include a hardening term and a recovery term in the 
work hardening state variable evolution equations. This form has been 
successful in modeling the Inelastic response of many isotropic metals. 
In this formulation an adequate data base for evaluating a recovery term 
did not exist and the term was not included. Also at the temperature 
studied, recovery is probably not as important as at higher 
temperatures. Future work should include the Investigation of a thermal 
recovery term in the state variable evolution equations. 

6.1.4 Latent Hardening 

The work hardening state variable evolution Equations (2.15) and 
(2.21) include "self hardening" only. That is, the work hardening on a 
slip system is assumed to be effected only by slip on that system. The 
hardening of inactive slip systems by intersection with active slip 
systems or "latent hardening" is generally considered to be an important 
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part of the theoretical basis for hardening in single crystal 
plasticity. References [64]-[673. In a recent review, Reference [773, 
the ratios of latent hardening to self hardening for single crystal 
aluminum and copper are reported in the range of 1 to 2. 

There are two methods for determining the effects of latent 
hardening mentioned in the literature. The first is to plastically 
deform a large crystal oriented for single slip, cut the crystal into 
specimens oriented to activate previously latent systems and compare the 
response of the initially latent systems to that in the original test. 
The second method involves measuring the crystal axis rotations in a 
tension test. During finite straining the slip direction rotates toward 
the loading direction and eventually a second slip system, or 
"conjugate" system, becomes equally stressed. If the latent hardening 
of the conjugate system is identical to the self hardening of the 
primary system it will slip equally under the same stress. If the 
latent hardening on the conjugate system is greater than the self 
hardening on the primary system equal inelastic straining will not be 
achieved until the stress on the conjugate system is greater than that 
on the primary system. 

Since Rene N4 i3 a high temperature alloy and the specimens are 
coated to prevent environmental degradation, the testing procedures 
outlined above would be further complicated. It seems that a viable 
alternative testing procedure for investigating latent hardening would 
be to run a combination of proportional and nonproportional tension- 
torsion tests. 
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For small strains and proportional loading, latent hardening will 
have little or no effect on the material’s response or the predictive 
capability of the constitutive model. For nonproportional loading or 
large strains the effect of latent hardening for Rene N4 is probably 
important but cannot be investigated with the current data base. 
Extension of the present theory to large strains and/or nonproportional 
loading will require consideration of latent hardening effects. 

6.1.5 Cyclic Hardening 

Nickel base single crystal alloys experience cyclic hardening 
resulting in up to a 1 0J increase in stress between initial and 
stabilized loops in a strain range controlled test, see Figure 6. 
Because of the limited data base this behavior has not been modeled. 
One way to model this effect is to include an additional term in the 
state variable hardness equation. A commonly used measure correlated 
with cyclic hardening is accumulated inelastic work which could be the 
independent variable in the cyclic hardening term. 

6.1.6 Other Response Characteristics 

There are a number of other response characteristics not modeled 
in this study. For example, slip bursts are observed at some 
orientations but have not been considered analytically. This effect is 
of little importance and in fact would present serious numerical 
problems if these effects were included in the constitutive model. 

Anelastic recovery, or negative creep strain rates at a positive 
stress following a stress drop, are commonly observed in nickel base 
superalloys. For Isotropic materials this behavior is usually modeled 
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by including a back stress tensor. A similar approach could be used for 
single crystal materials by including a back stress term on each slip 
system. In addition, the back stress model is important for isotropic 
modeling because it allows the direction of the deviatoric stress and 
inelastic strain rate vectors to be different and to vary as a function 
of the deformation. This effect has been observed and documented in the 
literature for a variety of multiaxial loading conditions. The current 
formulation does not force the stress and inelastic strain rate vectors 
to be parallel; however, there are no multiaxial tests of nickel base 
single crystal alloys to evaluate this effect. 

Tertiary creep is not predicted by the current model. The usual 
technique for Isotropic materials to account for tertiary creep is by 
Including a damage measure in the constitutive model. A similar 
approach could be used for Rene N4 by including another state variable 
and evolution equation. However, since tertiary creep occurs at strains 
above 2 % to 3%, there is little need to model tertiary creep for turbine 
blade and vane applications. 

6.2 Summary of Data Base Requirements 

For the nonlinear finite element code with the Rene N4 
constitutive model to be a useful tool it is necessary to develop the 
data base to generate the required material parameters. In addition, it 
is necessary to evaluate the model at other temperatures . The 
development of the data base can also be used for refining the 
constitutive model. In fact, part of the value of the constitutive 
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model is to provide a framework for alloy development and for 
understanding experimental results. 

The temperature environment for turbine blades and vanes varies 

from ambient up to about 1200°C. Inelastic behavior during operation 

would occur mostly at temperatures over 550°C. The test data for the 
inelastic response of the material should therefore span the range from 

about 550°C to about 1200°C. Considering the variation in response with 

temperature, it is desireable to obtain data at about 100°C intervals. 
A hypothetical test matrix at a single temperature is presented in Table 
6. The matrix is based on the tests required to determine the material 
parameters as discussed in Chapter III. The data base must be chosen to 
activate the octahedral and cube slip systems separately. Further , the 
tests should fully span the space of the sterographic triangle in both 
tension and compression. The [001] orientation is favored (half of the 
test3 are in this orientation ) since this direction is coincident with 
the radial direction of the blade or vane. In addition, tests should be 
run at other orientations to verify the model. The matrix in Table 6 is 
only for the model as presented in this report. The testing should also 
include a number of fatigue tests to evaluate cyclic hardening and 
further verify the model. It is not expected that the full test matrix 
should be run at all temperatures; however, the full matrix should be 
run at least one temperature above and below the critical temperature, 

about 800°C, since the deformation mechanisms appear to be different. 
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TABLE 6 ~ CONSTITUTIVE MODEL TEST MATRIX 
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TESTS IN THE [111] ORIENTATION ARE FOR EVALUATION OF CUBE SLIP CONSTANTS, TESTS 
IN THE OTHER ORIENTATIONS ARE FOR EVALUATION ON OCTAHEDRAL SLIP CONSTANTS. 



6.3 Conclusion 


The constitutive model developed in this study successfully 
predicts an Important part of the observed material behavior. A method 
has been developed for the derivation of material constants and no 
special or difficult tests are required. 

The implementation of the constitutive model in a nonlinear finite 
element code has been accomplished. It is computationally more 
intensive than comparable rate dependent isotropic models or rate 
independent anisotropic models, but it is not significantly more 
expensive. Some of the features in the finite element code enhance its 
usability and decrease its cost. 

The goal of developing a practical design tool for Rene N4 gas 
turbine engine components has substantially been achieved. 
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APPENDIX A 


EICSNSPACE APPROACH 

In the early stages of this research a considerable effort was made 
to further develop and apply the unified constitutive theory proposed by 
Stouffer and Bodner in 1979 C 38 3 • The theory as proposed included 
initial or deformation induced anisotropy. 

The form of the deformation rate equation is an extension of the 
Prandtl-Reuss flow law to anisotropic inelasticity. The rate of 
inelastic deformation tensor is related to the deviatoric stress tensor 
by a fourth order linear transformation whose components are functions 
of the stress and state variables. The constitutive equations are an 
extension of earlier work by Bodner and Partom [39], [40] and do not 
employ the use of a yield criteria or separate loading and unloading 
conditions. The state variables are introduced to characterize the work 
hardening of the material due to inelastic deformation. These variables 
are shown to transform as a fourth order tensor and are a central part 
of the anisotropic model. 

The tensorial structure of the constitutive theory is typical of 
most anisotropic constitutive models, namely 

e lJ “ A ijkl°ki (A,1) 

where is the strain or strain rate, is the Cauchy or deviatoric 


stress and is a fourth order tensor that characterizes all of the 

material properties. The components of ^ nay be constants as in 


linear elasticity or may be functions of stress and the previous history 
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as required by the particular choice of material and deformation 
mechanisms. Frequently Equation (1.1) is written in a six dimensional 
vector space as 


E a ” *a8°8 


(A. 2) 


Where a, 8 - 1,2,... 6 and and are written as vectors. In 
Reference [38] these vectors were defined as — a i i * °2 “ °22 °3 ” 


° 33 * - y/ 2~ <T 23’ °5 


/2~o^, ffg - /2~ o^ 2 and similarly for the strain 


tensor e^. 

The coupling between stress and deformation is obvious in Equation 
(A. 2). That is a single stress, 0 1 for example, will produce six 


components of deformation. One approach to modeling is to transform 
Equation (A. 2) into diagonal matrix; that is 

» 


* * 

X a 
act a 


(A. 3) 


* • » 

with no sum on a, e and 0 are the transformed variables and X are 

a a aa 

the eigenvalues of X .. This approach has the advantage that the stress 

QD 

and deformation are completely uncoupled; that is, one component of 
stress produces only one component of deformation. This could be 
particularly convenient in computational and experimental work. 
However, the major complication is that the eigenvectors are expected to 
be functions of X and, in general depend on the deformation history. 



The philosophy behind the proposed approach is to specify the 

eigenvalues and eigenvectors rather than the components of X a directly. 

as 

In Ref. [38], it was assumed that there existed a class of materials for 

which the eigenvectors are Invariant. 

Let us begin by calculating the eigenvalues and eigenvectors for an 

inelastic compressible isotropic material with a constitutive equation 

as described by Equation (A. 2). Designate e a3 the inelastic strain 

a 

rate and a as the Cauchy stress. The tensor X „ is assumed to be 
a aB 

symmetric, and since the components are real, then X^g will have real 

eigenvalues. Rice [41], Ponter and Leckie [42], Ponter [43] and 
Abuelfoutouh [44] have shown that the strain rate can be expressed as 
the derivative of a potential function with respect to the stress for 
history dependent polycrystalline materials. This result can be used to 
establish symmetry in X fl . The constitutive equation for an isotropic 

material in a six dimensional vector space can be written as 
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(A. 4) 
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The eigenvalues 


using two independent material functions A^ and A^. 
are 
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(A. 5) 


* * 

where e and o are the transformed strain rate and stress vectors in 
a a 


* t 

the eigenvector basis. The eigen configuration is defined as A • Q AC 


where Q is the proper orthogonal transformation that can be written as 

the matrix of the components of the eigenvectors. For this case the 
eigenvalues are not all distinct and the six eigenvectors are not 
unique. The only unique eigenvector is in the one direction. The other 
five eigenvectors are arbitrary, however, they should form an 
orthonormal basis. A typical orthogonal transformation, written with 
eigenvectors arranged as the rows of the matrix, can be defined as 
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A ^3 S3 

- " - o o o o 

/2 ✓ 2 


9 


1 I ' £ 

vi Jk Jk 


0 


0 


(A. 6) 


0 0 0 
0 0 0 
0 0 0 


1 0 0 
0 1 0 
0 0 1 


Using Q, the components of the transformed straip rate and stress 
vectors can be calculated in terms of the untransformed components. The 


eigenvector (1//3, 1 //3 » 1 //3» 0, 0, 0) i3 the hydrostatic line in 
stress space and the other two normal stress eigenvectors are in the 
deviatoric stress plane. Budiansky and 0'Connel [45] recently developed 
a method to evaluate the properties of isotropic polycrystalline 
material using this transformation. 
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A. 1 Cubic Materials 


Anisotropic materials with cubic symmetry have three independent 
material parameters X^, X 12 and X^. The constitutive equation, [8], 

can be written as 




in the principal directions of the material. This reduces to isotropic 
response when a relationship exists between the normal and shear 
components; i.e. X^ ■ ” * 12 ' 

The analysis for the eigenvalues and eigenvectors is identical to 

* 

that for Isotropic materials. The eigenvalues X , are X^ 1 + ^12' 

*11 “ *12’ *11 ' *12* *i»i 4 ’ *H 4 » x 44 for 0-1 to 6, respectively. 
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A. 2 Transverse Isotropic Materials 


Since a specific application is to metal inelasticity, it is 
reasonable to Impose the condition of incompressibility on the material 
response. This requires that - 0 in Equation (A.1) for arbitrary 


values of stress history and current loading. This result was derived 

by Olzak and Urbanowski [46] and Hill [47] for the general case of 

plastic flow with an anisotropic yield criterion. Introducing the 

deviatoric stress components, s, . - a.. - (1/3) tr 6,., into Equation 

lj lj mm lj 

(A.1) and imposing the incompressibility and symmetry conditions 


X lik£ “ X k£ii ’ 0 


(A. 8) 


gives 


E iJ “ *ijk£ s k£ + X ijkp ( 1 / 3 )ff mm 6 k£ " X ijk£ S k£ 


(A.9) 


Thus, for incompressible Inelastic material response let e Q and s q 

denote the inelastic strain rate and deviatoric stress vectors 
respectively, in a six dimensional space. In this space, the 
restriction of incompressibility can be written as six scalar equations 
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(A. 10) 
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The incompressibility condition specified in Equation (A. 10) can be 
incorporated into the r epresentation for a transverse isotropic 
materials [48] to yield 



Equation (A.11) was written with the coordinate axes in the three- 
dimensional space oriented parallel to the principal axes of the 
orthotropy of the body. In this case, coupling does not exist between 
the shear and normal components of the compliance matrix * a g- Again 

applying an elementary analysis gives three unique eigenvalues in the 
form 
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The eigenvectors can be arranged to give the transformation specified by 
Equation (A. 6). The eigenvectors for isotropic materials were defined 
to coincide with the eigenvectors for incompressible transverse 
Isotropic materials. 

A. 3 Incompressible Orthotropic Materials 

The material matrix X for orthotropic materials with the condition 

of incompressibility can be written as 



where e and s are again the components of the inelastic strain rate 
a a 

and deviator ic stress vectors. 

Solution of Equation (A. 13) for the eigenvalues yields six distinct 
eigenvalues, hence the eigenvectors are all unique. The 

» 

incompressibility condition requires >• 0 and the first eigenvector e 1 

- ( 1 //3, 1 //3 , 1 //3, 0, 0, 0) as before. However, the second two 
eigenvalues cannot be directly calculated; i.e. 
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(A. 14) 


A 22 * x 33 “ (A 12 + A 23 + A 31 ) - [ ' A 12 + A 23 + A 31 * A 12 A 23 


A 23 A 31 " A 31 A 12- 1 


1/2 


This further implies that the second and third eigenvectors, e 2 and e^, 
are functions of the tensor components *23 anc * A 31 * Since ttle 

matrix can be partitioned and analyzed as two independent three space 

A A A 

problems, the eigenvectors e^ , e 2> and e^ are all independent and can be 

A 

arranged as an orthonormal basis. Since e^is fixed for all 
incompressible orthotropic materials under any deformation history, the 

A A A A 

position of the e 2 and e^ can be defined by a rotation of e 2 and e^ 

A A A 

about the e 1 vector. Let 0 denote the position of e 2 and e^ relative to 

the vectors ( 1 // 2 , - 1 // 2 , 0 ) and ( 1 //?, l/£ - 2//S) in the three 
dimensional subspace, respectively. The first three components of the 
three eigenvectors can be written as 
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'1 


*S e - 
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1 
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0 - 
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0 

1//3 

1//3 

1//3 

cos 0 

sin 9 

\/2 

-1//2 

0 

-sin 9 

cos 0 

1//6 

iMT 

-2/S6 

J 




(A. 15) 
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The vectors e^ and e^ are still In the devlatorlc plane In the subspace. 

A A 

The angle 0 gives the history dependent position of and e^ and the 

representation reduces to transverse isotropy when 0 » 0. The 
transformation matrix Q, in Equation (A. 6), becomes a function of 0 for 

incompressible orthotropic materials. 

The angle 0 can be evaluated in terms of the components of X in the 

T * 

initial configuration. Noting that X - QXQ and that X^- 0, the 
calculation of X. 1 in terms of 0 gives 


tan 20 


6 , *12 " *23 . 

VTI 2A 12 - *13 ~ X 2 3 


(A. 16) 


6 r *22 ”*11 , 

“ IpT r— r — > 

/?2 33 11 22 

A. 4 Conditions Necessary for Stationary Eigenvectors 
A. 4.1 Transverse Isotropic Materials 

The minimum condition for constant eigenvectors can easily be 
established. If the transformation matrix §, defined by Equation (A. 6), 


is used to transform the matrix, X, 
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(A. 1 7) 



for an arbitrary compressible transversely isotropic material, the upper 

T 

left quadrant of the transformed matrix [QXQ ], becomes 


3 ( 2 A 11 + 2 A 12'” 4 A 13 +X 33 ) ° 


(X «-X -X ) 

/18 11 12 13 33 


A 1 1 -A 1 2 


^ (A 11 +A 12 _A 13" A 33 ) 


3 (A 11 +A 12 “ 4A 13 +2A 33 ) 


(A. 18) 


# 

This implies that the off diagonal term X^, must vanish if Equation 
(A. 18) is in the eigenspace. Thus, the condition 
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♦ X 


12 



(A. 1 9) 


is necessary to produce the constant eigenvector result. 

A. 4.2 Orthotropic Materials 

The compliance matrix for compressible orthotropic materials has 
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nine independent material parameters (X^, i 22 * * 33 * * 44 * * 55 . * 55 . ^ 2 * 

* 

*23’ *31 ^ 3ln 8 9 in Equation (A. 6 ) to transform A into the space A - 

7 * # * » 

QAQ , also implies that X 12> A^ and ^ must vanish if X is In an 

eigenspace. Thus 


* <2 ’ “it * V * “13 - V ’ ° 
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(A. 20) 


/T5 x 


31 


- X 


11 


*22 ’ 2 *33 + 2 *12 " *23 “ *13 “ ° 


are the necessary conditions for the eigenvectors to remain stationary. 
The first two equations of (A. 20) can only be satisfied if each of the 
quantities in brackets are zero. Thus 

*11 " *22 

(A. 21) 


*13 * *23 


and the third equation becomes 


*11 ' *33 * *12 *13 * ° 


(A. 22) 


These conditions are the same as that for transverse isotropy as 
expressed in Equations (A. 17) and (A. 19). Hence, materials with 
symmetry properties defined as orthotropic can never satisfy the 
stationary eigenvector condition. 

The consequences of this result can be further explained by 
considering an example calculation with the restriction that the 
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eigenvectors are fixed (this is the case presented in Reference [38]). 
Let us examine the Inelastic response of a strain hardening transverse 
isotropic metal. Suppose the hardening is calculated by modifying the 
eigenvalues in Equation (A. 12), and let e designate the inelastic strain 

I * * 

rates, e . Let X 22 and X^ be the eigenvalues after an arbitrary 

deformation and transform the eigenspace into the physical space 
according to Equation (A. 6 ). 

The resulting flow equation for the normal strain rates is given by 
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> (A. 23) 




Thus, the transverse isotropy assumption is maintained for any choice of 

* * 

X ^2 an< * *33 provided the transf ormation 0 remains constant. As a 


consequence, the material hardening computed in the eigenspace will also 
retain the transverse isotropy assumption. This result is 
unacceptable for general load histories 

The eigenvalues and eigenvectors must be functions of the material 
hardening for general load histories. Relating these functions to 
actual material response becomes a difficult, if not impossible task. 
In fact, the mathematical structure bears little relationship to the 
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physical metallurgy of the problem and merely confuses the situation 
rather than leading to any significant simplification. 
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APPENDIX B 


DEVELOPMENT OF THE KINEMATIC EQUATIONS 
This Appendix contains explicit representations for the local and 
global stress and inelastic strain rate components. The sign convention 
for positive slip directions on the (111) octahedral planes are shown in 
Figure B1 . The positive unit normals to the four octahedral planes are 
given by 

n 1 - 1//3 (I + I ♦ k) 

n 2 - 1//3 (- I ♦ I - k) 

n 3 - 1//3 (i - J “ i<) , and 

n 4 — 1 V 3 (- I - J ♦ k) (B. 1 ) 

where I, J, and k are the unit vectors In the principal material 

directions. The positive sign convention for the [1 0 1] directions on 
the octahedral planes are: 

1 11 - 1/V2 (I 

1 12 - w/2 (- 
5 13 - 1//2 (I 

1 21 - 1//2 (I 

1 22 - \/Jz (I 

1 23 - 1//2 (j 


- k) 

j + k) 

- J) 

- k) 

♦ j) 

♦ k) 
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PLANE 1 


PLANE 2 



PLANE 3 PLANE 4 



FIGURE Bl. OCTAHEDRAL PLANES AND SLIP DIRECTIONS 
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1//2 (I ♦ J) 



1 32 - l/vG (- J ♦ k) 

1 33 - 1/72 (I ♦ k) 

1 41 - 1/V? (j ♦ i<) 

1 42 - 1//2 (I ♦ it) . and 

1 43 - 1//2 (I - J) (B. 2) 

Substituting Equations (B.1) and (B.2) into (2.6) gives the inelastic 

strain rate tensor due to shearing in the [1 0 7] directions on the 
octahedral planes as: 
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Ill 



The Inelastic strain rate vector is relative to the principal material 
axes. Substituting Equations (B.1) and (B.2) into (2.9) gives the shear 

stress in each of the [1 0 i"] directions on the octahedral planes as: 



The stresses are relative to the material axes. 
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The positive sign convention for the Cl 2 1] directions on the 
octahedral planes are: 

I 11 - 1//S (- I ♦ 2J - ic) 

A 12 - \/JZ (2l - J - k) 

A 13 - 1//5 (- I - J ♦ 2k) 

I 21 - 1//5 (i + 2J ♦ k) 

A 22 - 1//S - (I - J - 2k) 

I 23 - l//£ (- 2l - J ♦ k) 

1 31 - 1/^? (- I ♦ J - 2k) 

1 32 - (21 ♦ J * k) 

A 33 - \/JZ (- I - 2j ♦ k) 

A 41 - 1 //? (- 2i ♦ j - k) , 

A 42 - 1/^ (I - 2J - k) , and 

A 43 - 1//6 (I ♦ J ♦ 2k) (B.5) 

Substituting Equations (B.1) and (B.5) into (2.9) gives the shear stress 

in each of the [T 2 7] directions on the octahedral planes as: 
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tension/corapression asymmetry for octahedral slip through the core width 
effect. 

The sign convention for positive slip directions on the (001) 
cube planes are shown in Figure B2. The positive unit normals are given 
by: 

-1 7 

n ■ i , 


-2 7 

n - J 


and 
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PLANE 1 



PLANE 2 



PLANE 3 



FIGURE B2. CUBE PLANES AND SLIP DIRECTIONS 
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n 3 - k (B.7) 

The sign convention for positive slip directions on the cube planes 

are: 

1 11 - )/S2 (] ♦ iO 

1 12 - \/Si (j - k) 

1 21 - 1/^2 (I ♦ io 

1 22 - 1/72 (I - k) , 

5 31 - 1/72 (I ♦ J) , and 

I 32 - lA'? (- I ♦ J) . (B.8) 

Substituting Equations (B.7) and (B.8) into (2.1) gives the 
Inelastic strain rate tensor due to shearing on the cube planes as: 
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Substituting Equations (B.7) and (B.8) Into (2.9) gives the shear stress 
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The total inelastic strain rate tensor is 
contribution from each of the active slip systems 
Equations (B.3), (B.6), and (B.9). 


the sum of the 
as computed in 
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